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Abstract 

An account is given of how the 'box integrals', as used for one- loop cal- 
culations in massless field theory, appear in momentum-twistor geometry. 
Particular attention is paid to the role of compact contour integration in rep- 
resenting the Feynman propagator in twistor space. An explicit calculation 
of all the box integrals, using only elementary methods, is included. 



1 Introduction 

The introduction of momentum-twistor space in (Hodges 2009) has been greatly 
developed by Mason and Skinner (2009) and effectively incorporated into the rapid 
development by Arkani-Hamed, Cachazo and Cheung (2009) of the powerful Grass- 
mannian formalism for describing scattering amplitudes in supersymmetric gauge 
theory. This note supplies some detail concerning the way that momentum-twistors 
can be used to represent the basic box-integrals appearing in the one-loop ampli- 
tudes. In particular, a twistor-geometric integral construction of thirty years ago 
takes on new life when re-interpreted in momentum-twistor space. We shall not 
be concerned here with reformulating the actual one-loop amplitudes which are 
found through the composition of these box-integrals. This, and the extension 
to higher-order loop integrals, is the subject of very active current work, notably 
by Arkani-Hamed, Cachazo and their collaborators. This note only discusses the 
integration which underlies these investigations. It is closely connected with com- 
plementary work by Lionel Mason and David Skinner. Their parallel publication, 
Mason and Skinner (2010), explores the connection of this integral formalism with 
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the AdS formalism and Wilson loops, and carries it further towards the evaluation 
of actual loop amplitudes. 

2 The box integrals 

The 'box integrals' that concern us are of the form 

f 1 d 4 p , 

J ((p - xx) 2 + ie)((p - x 2 ) 2 + ie)((p - x 3 ) 2 + ie)((p - x A ) 2 + ie) 

where the x±, x 2 , x 3 , X4 are 'region space' momentum parameters, characteristic of 
planar diagrams. These parameters express the conservation of momentum, by 
ki = Xi — Xi-i, and give rise to the concept of dual conformal symmetry. It is 
necessary to treat the non-degenerate case where all the fcj are non-null, and the 
various degenerate cases when some or all of them are null. 

The ie in the integral is the conventional indication that the real p-integral must 
be deformed into the complex in such a way as to avoid the poles of the integrand 
according to the Feynman prescription. But in what follows we shall actually 
discuss something slightly different and more general, namely 

/ ((p - Xl ) 2 - WKP ~ *2? - /i 2 )((p - x 3 ) 2 - WKP - - /i 2 ) ^ ' (1) 

where fx is a non-zero (complex-valued) parameter with the dimensions of mass. 
We shall be primarily interested in the limit when fx 2 — > from the correct direction 
in the complex plane, as this will recover the Feymnan prescription for massless 
field theory, but the actual integrals we evaluate will have a non-zero fx 2 parameter. 
Such integrals are finite. This notation also has the advantage of avoiding confusion 
with the quite different e used in dimensional regularization. 

There is a further generalization, in which each factor has a different fa mass 
parameter. Our basic geometrical setting permits this generalization, but the 
technical computation of integrals is more complicated. This question is briefly 
addressed in Section 8 below. Another reason for a focus on a single parameter 
parameter fx is that this should suffice to capture the information which is usually 
expressed with dimensional regularization methods. Our /^-dependent results may 
be regarded as transforms of the results as obtained in terms of a dimensional 
regularization parameter e. 

A further distinction must be drawn between the box integrals evaluated in this 
note, and the box functions generally tabulated. The box integrals A, as functions 
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of momenta, always take the form A = //A, where the numerator / is double- 
logarithmic and the denominator A is purely algebraic. There are good reasons 
for regarding the / as dimensionless multiplying factors, for seeing these factors 
as capturing the real content of the 1-loop integration, and for tabulating them 
as the 'box functions.' On the other hand, we shall exploit the fact that except 
in the most degenerate cases, / vanishes when A vanishes, and that the integral 
A itself is then finite. So we shall state the results in terms of the integral A, 
thereby including this special situation. The corresponding / can be obtained by 
a trivial removal of the rational denominator. But there is a subtlety: the box 
functions as conventionally defined absorb an important extra factor of | which 
is needed to make the correct connection between the scalar box integrals and 
the actual loop amplitudes. Thus the box functions F are actually of the form 



F=\f= |AA 

3 Wick rotation and contours in twistor space 



As is very well known, the Feynman contour prescription is equivalent to making a 
'Wick rotation' of the p-integration into the complex, so that the time- component 
of p runs along the imaginary axis. The integral thus becomes a Euclidean space 
integral. Adopting this point of view, we start by evaluating ([I]) in the special case 
of coincident x,. Without loss of generality, X, = for each i. The box integral 
becomes 



Let p° = it,p l = x,p 2 = y,p 3 = z, then take hyperspherical coordinates in M 4 for 
(t, x, y, z) and the integral can be evaluated as 



where the 2tt 2 factor comes from the 3- volume of a unit S 3 . The overall sign 
depends on contour orientation, equivalent to a choice of Feynman or anti- Feynman 
prescription, and in what follows we shall neglect it. 

The twistor translation of this idea turns out to give a striking example of twistor 
geometry. However, the first step only goes half-way to twistor space; we only go as 
far as the representation of complexified Minkowski space CM by CP 5 . Of course, 
this representation involves the compactification of CM, and its transformation 
under conformal transformations, which we shall address shortly. 




(2) 
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The underlying relations between the conformal group C(l, 3), the SO(2, 4) acting 
on CP 5 , and the SL(4, C) acting on twistor space, go back to the very origins 
of twistor theory in Roger Penrose's first work. The reader is referred also to 
the paper by Mason and Skinner (2009), which gives an introduction to twistor 
geometry using CP 5 as a bridge from Minkowski space. For present purposes, 
the most important feature of the correspondence is that the points in CM are 
represented by simpleskew bi-twistors Q a/3 in CP 5 , i.e. Q a/3 such that Q a @ = A^ a B^ 
for some twistors A a ,B^. This condition is equivalent to e. a ^Q a ^ QP & = 0, or to 
the rank of Q, considered as a linear transformation, being 2. Geometrically, 
this condition characterises the Klein quadric on CP 5 . Points of CP 5 not on the 
quadric correspond to complexified spheres on CM. These are essentially the very 
hypersurfaces of form (p — x) 2 — [i 2 = which appear as poles in our integral. As 
we shall see, the linear structure on CP 5 is highly advantageous in the handling of 
these poles. 

Ase aj3l s acts as an inner product structure, it is natural to introduce the notation 
X.Y for ^e a p 1 sX a ^Y 1& . We shall also need the special element I a/3 , corresponding 
to the vertex of the null cone at infinity, and the origin bi- twistor O a/3 . The infinity 
bi-twistor I a P is defined non-projectively, and knowing its scale corresponds to 
knowing the metric on CM. The relation is given by: 

= -'(fife- (3) 

We may let simple skew bi-twistors ,X% ,X^ correspond to the points 

xi,X2,X3,X4 in region space. We may also regard the variable p as an internal 
region and give it bi-twistor coordinates P a/3 . In the = 0, as studied above, 

the integrand of ^ then translates into 

(LP O.I) 4 



((20 + O.V/).P) 4 

The bi-twistor (20 a/3 + O.Ip, 2 I al3 ) is an example of a non-singular element of the 
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So far it would appear that twistors are being used only to express these bi- 
twistors; equivalently, only the lines of projective twistor space actually play a 
role, and not the points or planes. Moreover, since the box integral is a completely 
scalar structure, there might appear to be no reason for the spin-structure implicit 
in twistor space to be relevant. It is only when we seek to translate the differential 
form and the contour in (J2|, that we find twistor space coming naturally into the 
picture. 
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Suppose we do indeed resolve a bi-twistor into two representative twistors, and so 
write P af) = Z a V (3 -V a Z' 3 where Z a = (ip AB 'a B ,,a A/ ), V 13 = (ip BA ' t a , ,t b >). 

Then d 4 Z A d 4 V = i(aA'T A ') 2 d 4 p A d 2 a A d 2 r, where o~a>t a ' is the same thing as 
I a pZ a VP . (Note that here and in what follows we shall use non-pro jective twistor 
spaces and differential forms for the calculations, but this is only conventional; at 
each stage the integrals could be re-expressed in projective spaces if desired.) 



Now consider the eight-dimensional integral 

d 4 ZAd 4 V f id 4 P Ad 2 a Ad 2 r 



(2(O.I)-iO a p + Li 2 I a p)Z"Vfi) 4 J (p 2 - /i 2 )VA'r 



-A'\2 



(4) 



where the contour is defined by letting p run over the (non-compact) 1R 4 as defined 
above, and the spinor integral to be taken independently of p over an S 3 x S 1 . 
This may be done by taking an 'antipodal' or 'anti-diagonal' contour where ta> = 
t AA a a-, for some Hermitian non-singular t AA . It gives a factor (2iri) 3 . Thus the 
result (up to sign) is 



7T 
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r(27ri) 3 , 



6(/x 2 ) 2 

yielding a representation of the Feynman integral by a twistor integral. But it 
might be considered that genuine twistor structure has still not played any part. 
What matters in the integral is only the integration over the P a/3 . The spinor 
integral is a trivial factor; it merely integrates over the different ways in which 
pa/3 j g re p resen ted as a skew product of a Z a and a V 13 , and these representations 
do not have any significance. We seem to have added four extra dimensions, and 
then to have integrated them out, for no purpose. 

Twistor-geometric structure emerges only when we compactify this integral. To 
achieve this compactification, first let t AA be the unit vector in the O-direction. 
We can then define the Euclidean dual of a twistor Z a = (u A , tta 1 ) to be the dual 
twistor Z a = (t AA 'iiA,tAA'ti A ')- 

We can also define the Euclidean norm of a twistor Z a = (u A , ita') to be given by 
\Z\ = VZ a Z a = VoJ A tu A 't A A> + Tr A in A t AA '. = Vl^T + l^ 1 ! 2 + M 2 + M 2 . 
Next we consider the space of normalised twistors described by: 

(ip AB 'a B ',o-A') 



\(ip AB 'a B >,a A ') 



where p runs over the same space as before. These satisfy \Z\ 2 = \uj°\ 2 + \oo 1 \ 2 + 
1 7r ' | 2 + ki'l 2 = 1, but do not fill out a complete S* 7 , because the twistors of form 
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(p A , 0) are absent. These, of course, belong to the null cone at infinity in CM. 
Now note that provided p ^ 0,p ^ 0, we also have 



\(p A ,ipAA>p A )\ 

where p AA ' = p AA ' /p 2 , and p A = ip AB ' o B '- So we may identify p ^ 0,p ^ as 
defining two coordinate patches, together covering an extended contour. The point 
p = gives a one-point compactification of the original Euclidean IR 4 , turning it 
into an S" 4 . Moreover, the addition of the twistors of form (p A , 0) means that the 
seven-dimensional space of normalised Z now fills out the complete S 7 . For every 
projective twistor in CP 3 , there is an S 1 of points in this S 7 , thus realizing the 
Hopf fibration. 

Note that the Euclidean dual of Z a is 

(t AA 'a A , -it AA ,p A ' B o B ) 



% 

\(ip AB 'a B ',cr A/ 



but that for p of our form, p coincides with the time-reversed p, which means that 
pA'B _ p ; 



^a'b _ „a'b _ 2 (p.t)t A ' B = t AA 'p AB ,t BB ' . Hence 



If we further define 
then we have 



ya 



= {t AA 'a A ,-ip*t B ,a B ) 
\(ip AB 'cr B ',cr A/ )\ 

t^AB' + B zz. +AA' z=. \ 
[ip t B ,(TB, t a A ) 

\(ip AB 'a B >,(T A >)\ 



The pair (Z, e l6 V) then gives an eight-dimensional set with S 7 x S 1 topology. Inte- 
grating over the coordinate patch p ^ amounts to exactly the same as @. The 
remaining points, with p — 0, are only a set of measure zero and make no differ- 
ence to the value of the integral. Thus, integration over this complete compactified 
contour is equivalent to the Feynman integration. But the topology is completely 
changed by the compactification: in the compact contour, the spinor integral does 
not factorise trivially, and the integration is truly twistor- geometric. 

Having defined this compactification of the contour, it is much simpler to work in 
twistor space without any reference to integration in Minkowski space. Indeed this 
was the original construction in (Hodges 1977), as described in the historical note 
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in section 10 below. The starting-point then was the fundamental twistor integral, 
as described for instance in (Penrose and McCallum 1972): 



Qd 4 Z Ad 4 W _ 
(2vrz) 5 J (Z a W a ) 4 ' u 

where the contour can be represented by a S 7 x S 1 . By the coordinate change 
W a = QapV 13 , we have 

6 d 4 Z A d 4 V 1 



(2™) 5 / {Q aP Z«V?y detQ' 

where Q is any non-singular linear transformation. If Q is antisymmetric, the 
determinant simplifies to the square of the Pfaffian and we have 



1 f6d 4 ZAd 4 V 16 



(2tu)5 J (Q a pZ«V?) 4 {e^ 5 Q aP Q 1& Y ' 
In particular, for Q a/S = 2(0.iy 1 O a ^ + /i 2 J a/3 , we obtain 

6d 4 ZAd 4 K 16 



(6) 



{2mf J ((2(0./)- 1 0^ + ^I af} )Z*VP) 4 {8fj. 



2\2 



(7) 



in agreement with (|4]). This is the fundamental compact twistor integral on which 
everything that follows rests. 

An important feature of the compactification is that the non-singularity of the 
form on the compactified contour is equivalent to the absence of an ultra-violet 
divergence. 

When the Feynman integration over momentum p is recast as this compact twistor 
integration of a holomorphic form, it becomes free from its original definition in 
terms of spacetime. The Feynman propagator is thus given a new geometrical 
interpretation, which is essentially defined by the CP 3 of projective twistor space. 
It is a striking fact that from this point of view, a loop integral has the shape 
of twistor space. A further freedom arises from the fact that the twistor-space 
contour, more correctly regarded as a homology class, need not possess any imme- 
diately obvious connection with the topology of an S 7 . In section 10, an example 
is given of this freedom of representation. 

The underlying structure is the quaternionic fibration of S 7 over S 4 , analogous to 
the Hopf fibration. That is, we consider an S* 7 as the space of pairs of quater- 
nions (gi,^) with |gi| 2 + |g 2 | 2 = 1- Then define equivalence classes on S 7 by 
Q2) ~ {q[, q'2) ^ 1i = IQiy I2 = Q.Q.2 f° r some unit-normed q. The quotient space 
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is homeomorphic to S 4 and each fibre is an S 3 . This fibration rests on the fact that 
the quaternionic norm satisfies \qr\ = \q\\r\, which itself can be expressed simply 
as the four-square identity for real numbersF] 



4 Generalised Feynman parameters 

From now on we can abandon the original p-space integration, and derive every- 
thing from the compact twistor contour integral: 

1 [Qd 4 ZAd 4 V 16 

( 8 ) 



This is valid for all non-singular anti-symmetric Q in CP 5 , i.e. those which are 
not on the quadric &*^ s Q a pQj6 — which corresponds to complexified Minkowski 
space. The difficulty in our programme, of course, is that we are interested precisely 
in the limit fi 2 — > where the limiting Q is on this quadric, so that this compact 
twistor contour no longer exists. 

We now address the original integral, with its four different Xj. Rather than find 
a new contour in twistor space, we think of moving the Xi apart from coincidence, 
while keeping the twistor contour the same. There are very well known techniques 
for doing just this, using Dirichlet averaging, in the form used extensively by 
Feynman. The principle is that of embedding the integral in a larger space and 
then exchanging the order of integration. Feynman's classic method uses the 
identity expressed in a symmetrical form by 

6 5(1 - «i - a 2 - tt3 - «4) , j i , 

daida 2 da 3 da 4 




Q1Q2Q3Q4 Jo Jo Jo Jo (ttiQi + CH2Q2 + a 3 Q 3 + a 4 Q 4 



which is an integration over a tetrahedron in M 3 . But we are at liberty to prefer 
the identity: 



6 da d/3 dj 

(9) 



Q1Q2Q3Q4 Jo Jo Jo (Qi + + PQ3 + 7Q4) 4 ' 

Here and throughout, integration to 00 is to be interpreted as a compact integral 
in a suitable (complex) projective space; so in this formula, a,/?, 7 are actually 



1 I am indebted to Lionel Mason for pointing out that the S 7 x S 1 contour must be connected 
with the Euclidean space by this fibration. Earlier work had used the S x S , noting that it 
gave the correct answer, without seeing the underlying reason for this correctness. 
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CP 1 parameters. The topology of the space is that of a cube in (CP 1 ) 3 . This may 
be considered to correspond with the Feynman tetrahedron by 

a = «i(l — ol\ — a>2 — as) -1 , ft = a2(l — oii — 0:2 — a^) -1 , 7 = 03(1 — cti — 0*2 — 0^3) 1 , 

a mapping which blows up the = 1 face into the three faces of a cube given 
by a = 00, (3 = 00, 7 = 00. Clearly, this parameter space breaks the symmetry 
which is respected by the Feynman tetrahedron, but as a direct product of three 
line intervals it may be more convenient for computational purposes. We shall also 
use a 3-volume with the shape of a triangular prism, and the identity 

1 _ [°° dx f 6 ( X + x ) 2 ( 1 - "i - "2) daida 2 

Q1Q2Q3Q4 ~Jo X Ja ((1 + x)(axQ 2 + a 2 Q 3 ) + (1 - an - a 2 )(xQ l + Q 4 )) 4 ' 

(10) 

where A is the triangular region < ai,0 < a 2 ,ai + a 2 < 1. This can be 
regarded as arising from taking the Feynman tetrahedron and replacing by 
x = 03(1 — «i — «2 — c^)" 1 as a parameter. This blows up one of the edges of the 
tetrahedron into a rectangular face. 

Such transformations of the Feynman parameters are nothing new, and indeed no 
new results will follow from the analysis in this note. The most general integral of 
this form, with different /ij, was calculated exactly by t'Hooft and Veltman (1979) 
long ago. We are merely reviewing the derivation of the main results in the new 
context provided by the compact twistor-space contour and the use of the linear 
structure of CP 5 instead of that of 1R 4 . 

The Feynman parameter method removes the need to do any integration over the 
original space (here replaced by a twistor space), once the fundamental integral 
has been done. If both the Feynman parameter space and the twistor contour 
are compact, then the change of order of integration can be rigorously justified. 
However, this depends upon the same twistor contour being used for each param- 
eter point, and this is harder to establish. The fact that the twistor contour is 
compact gives a very useful starting-point. Having established a specific contour 
for ([8]), for some specific Qo, the same contour must be valid for Q in some open 
neighbourhood of Qo, and in particular for a polyhedron in that neighbourhood. 
Restricted to this region, therefore, the method is rigorous. In practice we extend 
to all non-singular Q by an argument from analytic continuation, thus assuming 
in effect that the contour can always be moved around as required in the twistor 
space. Indeed we can go further and by moving Q parameters in a loop, deduce the 
existence of period contours. Whilst more careful mathematical study would be 
desirable, we are not seriously concerned about the validity of the results obtained 
by informal methods. 
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5 Evaluation of the 4- mass case 



The '4-mass' box integral is simply the non-degenerate case where all the external 
Ki are non-null, and so equivalently, all the Xj are non-null separated. In CP 5 
language, the corresponding elements satisfy Xi.Xj ^ 0. The corresponding lines 
in twistor space are all skew. We now wish to evaluate the box integral ([!]) in this 
case, and in the limit as /x 2 — > 0. 

We have established that by integrating over the direct product of a compact 
S 7 x S 1 contour in the twistor space, and a cube-shaped region in the Feynman 
parameter space, 

6d 4 Z Ad 4 V 



(2tu) 5 J (Q la pZ a VP)(Q 2a pZ a Vf ) )(Q 3a f ) Z a Vf } )(Q 4a f,Z a Vf ) ) 

dad/3d 7 



JO JO 



det(Qi + aQ 2 + /3Q 3 + 7Q4) ' 

So we may apply this to the case where = 2{Xi.I)^ 1 X^ + fi 2 I alS , and the 
Xi correspond to four skew lines in twistor space. Using the abbreviation xfj for 
-2(X i .X j )(X i ./X i ./)- 1 , the result is 

da d/3 d7 



OO POO POO 



Jo Jo (ax 2 12 + (3x 2 13 + ^x 2 4 + a(3xl 3 + ajxl 4 + (3jx 2 u - fi 2 (l + a + P + 7) 2 ) 2 ' 
where the connection with the standard kinematic parameters is given by 

2 ± 2 2 7.-" 2 2 _ 7. ' 2 2 _ 7>"2 2 j, r 2 

x 13 — Z i X 24 ~ s > X 14 ~ "l> X 12 ~ "25 • r 23 ~~ "35 X 34 ~ -*M ■ 

Here we take a short cut by taking the limit fi 2 — ¥ inside the integral. In doing 
so, we are neglecting to take proper account of how the limit fi 2 — > encodes the 
Feynman prescription. In effect, we shall find the amplitude function as a many- 
branched complex function, leaving to later the question of which branch is the 
correct one to take. 

The 7 integration is trivial and leaves 

da d/3 



OO POO 



'o Jo ( aX 12 + P X 13 + afiAl) { X \4 + «^24 + 

iarameters admit a 
which transforms the integral into 



The a and /3 parameters admit a simple rescaling u = ax\Jx\^ v = f3x\ 4 /x\ A , 



du dv 



Jo (^12^34 + VX 13 X 24 + UVX% 3 xl 4 ) (1 + U + V) 
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so that evaluation amounts to performing the double integral 

du dv 



oo roo 

t] ii or 

12 



o J {auv + bu + cv)(l + u + v) 

where a = x\^x\ 4 = K%Kf, b = x\ 2 x 2 u = K 2 K 2 , c = xl 3 x% 4 = st. Note that the 
ratios a/b,b/c, c/a are conformal invariants. 



Performing the integration over v leaves 

log ((ait + c)(l + u)/bu) 



au 2 + (a + c — b)u + c 



du. (13) 



We note that 



log((a M + c)/H ^ = r log((w + l)a/b) 



au 2 + (a + c — b)u + c J aw 2 + (a + c — b)w + c 



by w = c/au, so that the integral (13) is equal to 

21og(l + M ) + log (a / fc ) dM 
au 2 + (a + c — b)u + c 

In what follows we shall write (—ft) and (—ft) for the two roots of the quadratic 
au 2 + (a + c — b)u + c, so that nk = c/a and ft + k = (a + c — b)/a. We also write A 
for a(ft — ft), so that A 2 = (a + c — 6) 2 — 4ac = a 2 + b 2 + c 2 — 2ab — 26c — 2ca. 



In these terms the integral ( 14 ) becomes 

(21og(l + w)+log(a/6)), 1 J" (2dilog(l + u) + log(a/6) logw) , 

-du = - — - ® -. r- — du . 



a(u + k)(u + ft) 2iri J a(u + K)(u + n) 

by using the property of a contour surrounding the dilogarithmic cut (see equation 
(57) in the Appendix, where some basic properties of the dilogarithm are listed). 
Deforming the contour into circles around the two poles, we finally evaluate the 
integral as: 



2 dilog(l - re) - 2 dilog(l - re) + log(a/6) log(re/re) 

(15) 



A 

2dilog(l - re) - 2dilog(l - re) - log((l - re)(l - re))log(re/re) 



A 

This can be written in many other forms, because any permutation of {x\, X2, x$, x±) 
and hence of (a, 6, c) corresponds to a cross-ratio transformation generated by 
re — > 1 — re, re — > re" 1 . The dilogarithm is invariant, up to logarithmic terms, under 
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such cross-ratio transformations, and so numerous identities are available. Using 



one such dilogarithmic identity (55) it can immediately be written as 



2dilog(l - fiT 1 ) - 2dilog(l - - log((l - - k~ 1 )) log(«//e) 

A 

which is equivalent to the exchange of a and c. Averaging over these two expres- 
sions gives a longer formula which is equivalent to the expression preferred by Bern 
et al. (2004), in their equation (41). 

Another expression, longer but more elegant and manifestly invariant under per- 
mutations of (xx, X2, xs, X4), is given by 

1 J di\og(n/k) — dilog(/«//«) — dilog(l — k/1 — k) + dilog(l — k/1 — k) 
A \ +dilog(l - k" 1 /! - - dilog(l - k^/l - k' 1 ) 

(16) 



To establish this, write ( 14 ) as 



log((l + u)/(l - k)) + log((l + u)/{l - k) 

'/ w ^ du . 

a[u + k)[u + K) 

For the first logarithm, change the contour to run from to — K, where there is 
a removable singularity, and then from — k to 00, changing variables by a Mobius 
transformation to z = —(u + k)/(u + R). For the second logarithm, exchange the 
roles of k and k. The result is then immediate. This use of a Mobius transforma- 
tion based on a removable singularity is equivalent to proving Abel's fundamental 
identity for the dilogarithm. 

A special case arises if the roots are coincident, so A = 0. This condition, equiv- 
alent to y/a ± \fb ± \fc = 0, is closely analogous to Ptolemy's condition (that 
\AB\\CD\ ± \AC\\BD\ ± |AD||SC| = 0) for four points A,B,C,D in the Eu- 
clidean plane to be cocylic. Geometrically, A = implies that the four points 
lie on a complexified circle, or that that the four elements of CP 4 are linearly 
dependent. In this case the integral is still finite, and its value can be written 
symmetrically as 

log a log b log c 



From the uf-integral (12) one may read off the existence of three period contours, 



obtained by analytically continuing in a, b, c. These are given by 

2m / — -. du, 2ni / — r du , 

J au A + [a + c — b)u + c J bu z + [b + a — c)u + a 

f°° 1 

2-ki I — - — r -du, (17) 

J n cu 2 + (c + b - a)u + b K ' 
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which yield 



A A A 

again with special cases when A = 0. These three contours add to zero, modulo 
the double period which gives (27ri) 2 A _1 . It is the double period which corresponds 
to the 'leading singularity' of the amplitude, and also corresponds to the integral 
obtained by changing the Feynman propagator to a 5-function on the light-cone, 
i.e. putting the propagator 'on-shell'. This 'leading singularity' structure is very 
significant in the current development of loop amplitude theory. Analogous phe- 
nomena of period contours have appeared in twistor-theoretic literature since the 
very earliest days, though in a different physical context. A historical note is given 
in section 10. 

This period structure is relevant to a question which we have not yet addressed, 
namely that of how the Feynman prescription is translated into the correct choice of 
twistor contour, and thus into the correct branch of the dilogarithm. Here the work 
of Duplancic and Nizic (2002) emphasises a subtle breaking of conformal invariance 
in the amplitude. Although the dilogarithmic function, as a complete analytic 
function, is a function only of the conformally invariant ratios of a/b, h/c, the choice 
of which branch to take depends on knowing the individual variables s, t, Kf. Thus, 
the /x 2 parameter leaves behind a remnant of the infra-red divergence it regularises, 
even when the limit is finite. This is evident from the original definition of the 
integral, in which the Feynman prescription applies to the individual Xj. In the 
derivation of the dilogarithmic function we have neglected this prescription in two 
places: first when the limit fi 2 = is taken before doing the integral, and then 
in the rescalings u = ax\Jx\^v = fix\Jx\±. One may start with a configuration 
where all the Xi are spacelike-separated; in this case there is no ambiguity in 
the integral, nor in the rescalings, and u, v, w may all be taken to run along the 
real positive axis. But analytic continuation of the resulting amplitude to timelike- 
separated Xi requires knowledge of the Feynman prescription. Duplancic and Nizic 
express the choice of branch in terms of rules for adding on a logarithmic period 
function. 

It is worth noting that Duplancic and Nizic (2002) make use of the triangle integral 
in order to make their definition of the correct dilogarithmic branch. In our picture 
this triangle integral also has a natural geometric meaning: it simply corresponds 



to putting Q4 = / in (11 ) 
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6 Momentum-twistor parameters and 
transversals 



So far we have expressed the external parameters Xi in terms of elements Xi of 
CP 5 . But we can also separate these CP 5 elements into representative twistors. 
They can immediately be identified as the momentum-twistors as introduced in 
(Hodges 2009). 

One motivation for doing this is that when the box-integrals are combined with 
the helicity structure of actual gauge-fields, we shall certainly need such helicity- 
carrying twistor parameters. But even in the 4-mass scalar box-integral, where no 
helicity structure is apparent, the twistor-space structure is a useful adjunct. 

In (projective) twistor space, the four external Xi will correspond to four skew 
lines. The analysis of the 4-mass integral has shown that it behaves as a residue 
calculation, the poles being determined by k and k as the solutions of a (confor- 
mally invariant) quadratic equation. This structure reflects the existence of just 
two transversals to those four skew lines. The special corresponds to 

the coincidence of those two transversals. 




Figure 1: 4-mass transversals 

In Minkowski space, the transversals correspond to the solutions p ± of the equa- 
tions (p — X\) 2 — (p — x 2 ) 2 — (p — x 3 ) 2 — (p — X4) 2 = 0. In principle this is just a 
quadratic equation, but it is surprisingly difficult to write down a formula for p ± . 
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The reason is that giving a formula in Minkowski coordinates involves a fifth skew 
line in twistor space, namely the line corresponding to the point at infinity. This 
introduces many more algebraic invariants which are not actually relevant to the 
geometry of the transversal itself. 

The geometry is much simpler in CP 5 , where the four points Xj define a linear 
subspace. This has an orthogonal subspace (with respect to the inner product 
defined by € a p~,s), within which just two elements lie on the Klein quadric and so 
represent points rather than complexified spheres. 

In twistor space, the geometry is particularly elegant. Three lines, say those of 
Xi, X2-, £3, define a quadric, ruled by two families of lines, each line being transver- 
sal to every line of the other family. The line corresponding to 24 then meets 
the quadric in just two points. These two points define the two lines which are 
transversal to all four given lines. The construction is manifestly conformally 
invariant. 

Explicitly, the quadric is e^sex^X^ 1 X^ X X^ 'Z a Z a = 0; if P a and Q a are two 
points on the line x±, then there are just two (possibly coincident) roots z for 
the twistor Z a = P a + zQ a lying on this quadric; the transversals are then given 
by T±jip = Cap-yStXtipaXi" 'X% X Z±Z±. By choosing the basis twistors P a and Q a by 
reference to the line at infinity, one may obtain a formula for the transversal as a 
point in Minkowski space, but the properties of the transversals do not depend on 
knowledge of the infinity twistor. 

Each transversal is actually a CP 1 , and the values the cross-ratios of the 

points where the four lines intersect it. 

The transversal picture illuminates the structure of the degenerate cases where 
one or more of the external momenta is null. If k 2 = for some i, then the 
quadratic equation au 2 + (b — a — c)u + c has a root at 0, 1, or 00. Analytically, 



the dilogarithm in (16) is then divergent, and this reflects the fact that the fi 2 — > 
limit is no longer finite. In the next section we shall calculate the form of this limit. 
But the geometry of the transversals is still well-defined and very simple. 

Consider first the 3-mass case, with k\ = 0, which corresponds to X\ being null- 
separated from Xi, and so to the corresponding lines in twistor space having a 
common point. The other Ki are non-null, and so all the other lines in twistor 
space are still skew. There are still just two transversals to all four lines, but 
now they are distinguished in the twistor-space picture by a simple geometrical 
criterion: there is one transversal through the point common to X4 and x\, and 
one in the plane containing X4 and X\. 

This is quite different from the 4-mass case, where nothing in the geometry distin- 
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xi 




x 4 



x 2 x 3 



Figure 2: 3- mass transversals 



guishes one transversal from the other. This bifurcation into two types of term, 
governed by the two possible helicity representations, is very well known. When 
the box integral is combined with the gauge field structure in order to calculate 
actual amplitudes, summation over these qualitatively different terms plays an es- 
sential part. (Indeed this summation, including the important factor of 1/2, gave 
rise to the original BCF recursion relation.) Here we see a natural characterization 
of this bifurcation in terms of the geometry of transversals in twistor space. 

The momentum-twistor space parameters extend to the degenerate cases very sim- 
ply. As a matter of convention, we shall allocate names to them in such a way 
that Ki is associated with A twistors, K 2 with B twistors, and so on. In the four- 
mass case we may take X\, x 2 ,x 3 , x± to correspond to lines A 2 Bi, B 2 C\, C 2 Di, D 2 A\ 
respectively. (In applying these results to the actual calculation of n-point am- 
plitudes, we would in general have more than two twistors associated with each 
but allowing this generalization would not add anything to the analysis at this 
point.) 

For the degenerate cases, the momentum twistors associated with a null momen- 
tum simply merge into one. Taking the 3-mass case, we may consider the external 
momenta as parametrised by seven momentum twistors A, B\,B 2 , C\,C 2 , D\, D 2 , 
where X\, x 2 , x 3 , X4 correspond to lines ABi, B 2 Ci, C 2 Di, D 2 A respectively. Ex- 
plicitly, we then have 



t=(k 1 + K A ) 2 = ( Xl - x 3 ) 2 



2 (ABi) (CW 



(AB.CM 
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s = (h + K 2 y = (x 2 - x 4 y = -2 



o {AB 2 C 2 D 1 ) 



(B 2 C 2 ) (D 1 A) 



K -^-^- 2 u^Wrr etc - (18) 



where (PQRS) stands for e a p lS P a Q fi W S & and (PQ) for I a pP a Q ri . 

Although these statements have been made for twistor space, they could equally 
well have used the dual twistor representation. In the dual twistor picture, the 
geometry of the the transversals is interchanged, point becoming plane, and plane 
becoming point. 

One transversal, that passing through the point A, is given by 

T a (3 = ^X^Tr[a^l3]paTA B%C 1 A P D 1 C 2 . 

The other, lying in the plane D 2 ABi, is given by the dual formula. These formulas 
also apply to the more degenerate cases. In the '2-mass-easy' case, where C\ and 
C 2 merge into C, one transversal is AC and the other its dual. 

In the following calculations of the scalar box function, it makes no difference to 
the technical work, whether the external parameters are momentum twistors or 
dual momentum twistors. We shall continue by referring to a parametrization with 
momentum twistors, leaving the dual case implicit. 



7 The box integrals in the degenerate cases 

In the degenerate cases, our problem is that in the /i 2 — > limit, the bi-twistor Q 
is simple not just at the vertices of the Feynman tetrahedron, but on from one to 
four edges. These singular edges are the source of the logarithmic divergence as 
fi 2 -> 0. 

In principle, one could use this new geometrical setting to re-derive the results 
of t'Hooft and Veltman (1979) by integrating over a general tetrahedron in the 
parameter space. This would give the most general picture of the divergence due 
to the singular edges. But we shall restrict the analysis to the situation of equal 
ji, and moreover neglect 0(/i) terms, i.e. those which vanish as /i 2 — > 0. There 
are many possible representations of the results, but we shall organise them in a 
form which expresses the //-dependence entirely in powers of log(/x 2 ), and which 
facilitates comparison with standard expressions. 

We shall use entirely elementary methods (as opposed to Mellin transform meth- 
ods), with the intention that the location of suitable co-ordinates, and the various 
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devices involving splitting and recombination of terms, may be helpful in further 
elucidation of the underlying geometric concepts. 

We shall also give an exposition in terms of the standard kinematic variables 
s,t, Kf. This is purely for temporary convenience, in that it is easier when work- 
ing through the details of the integration in the different degenerate cases to be 
reminded by the notation of the physical setting of null and non-null momenta, 
and thus of which quantities have been set to zero before the /z 2 — > limit is taken. 
It also assists in making contact with standard expressions. But this notation is 
in a sense unfortunate and retrograde, because it conceals two important aspects 
of the theory given here. 

Firstly, each of these variables is defined as a twistor-geometrical invariant of the 



momentum twistors Z{ and the infinity bi-twistor J, by the relations (18), and 
the final value for the integral should likewise be considered as such a twistor- 
geometrical object, not as a function on standard momentum-space. In particular, 
the emergence of terms which do not involve J, and so are conformally invariant, 
is a vital feature of the theory. 

Secondly, the momentum twistors should all be considered on an equal footing, 
rather than some being associated with null and some with non-null momenta. 
In the application to the calculation of loop amplitudes, a summation is taken 
over all the ways in which the momentum twistors can clump together into four 
subsets, and any particular momentum twistor Zj thus plays many roles within the 
summation. The complementary paper (Mason and Skinner 2010) uses a notation 
which properly illustrates both of these features of the theory, but which is less 
well suited to the details of integration on which we now embark. 



7.1 3- mass case 

In what follows, we shall as a matter of convention take in all cases that k\ is null, 
i.e. (24 — Xi) 2 = 0, and hence that the corresponding (4-1) edge of the Feynman 
tetrahedron lies in the singular quadric. 

The 3-mass case is the least degenerate, with only that one edge of the tetrahedron 
on the singular quadric. In this case we can use the same cube-shaped Feynman 
parameter space defined by ([9]), and the integral (11) becomes 



Jo Jo 



dad/3d 7 (ig) 



{a(Kl + js) + P(t + jKl) + apKi - /x 2 (l + a + (3 + 7 ) 2 ) 2 ' 



Note that the \i only plays a role near a — (3 — 0, i.e. at the edge lying in the 
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singular quadric. Hence we can replace this integral by 

da d/3 d7 



COO /"OO 




JO JO 



{a{K\ + js) + (3{t + jK%) + apKl - /i 2 (l + 7 ) 2 ) 2 



o(a0 



= 7o (t + 7 ^ )( K| + 7S )-^i(i + 7 ) 2 d7 + 0( ^ (20) 

The /i 2 term in the denominator can be dropped. This is because the denominator 
could be expanded as a Taylor series in /i 2 , but because, for n > 1, (/x 2 ) n log(/i 2 ) 
itself vanishes with /i 2 , all terms above zeroth order may be absorbed into the 
O(fj). The effect is to leave: 

- \og{-^Kl/Klt) log(KjKifst) 



st - K\Kl 



I log(l + 7gj/Q + log(l + Is/Kl) - 2 log(l + 7 ) 
: 1 {t + 1 Kl){Kl + 1 s) d 7 + 0(M). (21) 



Using the dilogarithmic relation (57) again, the remaining integral in (21) can be 
written as 

1 I dilog(l + jKl/t) + dilog(l + 1S /K 2 2 ) - 2 dilog(l + 7 ) 
2™ / (H7«T^ d7 + 

dilog(l - KlKl/st)- dilog(l - st/fffif 2 2 )-2dilog(l - iff /s)+2 dilog(l - t/ATf) 
~ - K\K\ 



Collecting terms, using the dilogarithm identity (55), we obtain: 



A 3m (s,t,KlKln 2 ) 



1 



log(-/x 2 #f/ S *) \og{KlKl/st) + 2dilog(l - KjKl/sty 
-2dilog(l - Kl/s) - 2dilog(l - K\jt) ) . (22) 



st - K\K\ y _j hg 2 {K 2 /s) _ l log ^ K l/ t) 



This may be written more symmetrically as: 
1 



\og{-n 2 Kl/st) \og(K 2 Kl/st) + 2 dilog(l - K 2 K 2 /st) 

, , , v — ( dilog(l — Kl/s) — dilog(l — t/Kj)) }. (23) 

K * K 2 [ _ ( dilog(1 _ K y t) _ dilog(1 _ s/K 2 2)) 

This expression makes it more transparent that there is no pole when st = K%K%, 
the analogue of the condition A = in the four-mass case. 
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The full symmetry in (s, t) -h- (Kf, K^) is manifested by the longer expression: 

f \og{- ^KHstKlKl)-^) \og{KlKl/st) 
1 I + dilog(l - KlKl/ st) - dilog(l - st /KlKl) 
st - KlKl ] -( dilog(l - Kl/s) - dilog(l - t/KD) 
{ -( dilog(l - K 2 Jt) - dilog(l - s/KD) 

In particular, the integral allows ki not just to be null, but to be zero, so that 
X4 = x\. The integral reduces in this case to a triangle defined by the three 
distinct Xi, (The 7 integral is trivial, as s = Kl,t = Kf.) The calculation above 
may be considered as treating this configuration as a base and then varying away 
from it with a non-zero but null k±. 




7.2 2-mass-easy case 

Here k% is null as well as k%, so two opposite edges of the Feynman tetrahedron 
lie in the singular quadric. We now have only six momentum-twistors; i.e. C\ and 
C 2 have merged into C. There is considerable symmetry: not only is the result 
symmetric under (iTf,iff) -h- (s,t), but also under K\ K\, s -f-)- t. There 
are generalized Feynman parameters which do better at keeping this symmetry 
manifest, but the following method has the advantage that it requires no solutions 
of quadratic equations, and also that it extends to the more degenerate cases. It 



uses the triangular prism (10), under which the integral becomes 

dctidcWl — cxi — ^2) 



00 

dx 



J J a ((1 - "1 - a 2 )(<*i(s + K \x) + a 2 (Kl + tx)) - /i 2 (l + x)) 2 ' 

where the region A is the triangle {0 < a<i < 1, a\ + a 2 < 1}. For the next step 
we give new coordinates to this triangle by y = 2(ai + a 2 ),u = 4ai(l — a± — a 2 ). 
These have the effect of blowing down the edge 1 — a\ — a 2 = of the triangle 
into a single point. The result is to map the original tetrahedron into a figure 
which is a parabolic segment times a line interval, with the original opposite edges 
becoming two parallel edges at y — 0, y = 2. The integral becomes 

c r-2 ry{2-y) j M 

o dX l dy l Ws + Kjx) + (y(2 -y)- u)(K* + tx) - V(l + x)f ' (25) 

This is clearly symmetric in y -H- 2 — y, so we split [0,2] into [0,1] U [1,2], and 
change y to 2 — y in [1, 2]. The effect is to concentrate the regularization entirely 
into the corner at u = 0, y = of the integral 

,00 ,1 f y(2-y) du 

2 dx dy 



^0 ^0 



u(s + K\x) + (y(2 -y)- u){Kj + tx) - 4/i 2 (l + x)f 
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Now we use 2 = 2y + (2 — 2y) to separate this integral into two pieces: 

f°° f 1 ry( 2 -y) (j M 
2 Jo ^ Jo V ^ Jo Ha + K 2 x) + (y(2 - y) - u)(K 2 + tx) - V(l + x)) 2 W 

and 

00 r l ry( 2 -y) ^ u 



dx (2-2y)dy 
Jo Jo 



'u(s + K$x) + (y{2 -y)- u)(K 2 + tx) - V) 2 (l + x)) 2 ' 

(27) 



The first piece (26) simplifies because it is finite as [i — \ 0, and can be evaluated 
immediately as 

2 ^ , log( KjKj/st) 

lo g 4 — T/2 TS2 



st ta - KfifJ 



The second piece ( J27| ) simplifies because the numerator (2— 2y) is just the derivative 
of y{2 — y). Let w = y{2 — y) — u replace y, and it becomes: 



2 / da; 



o 



du dw 



A (u(s + Kfx) + w(ATf + tx) - 4/i 2 (l + x)) 



2 



where A is the triangle < w,u,w + u < 1. Elementary integration of (u, w) over 
the triangle yields 

r \og{s + K 2 x)/(s + K 2 x) - \og{K\ + tx)/(K 2 + tx) 
Jo (s + Kix)-(K! + tx) [ZJ) 



The expression ( 28 ) combines with ( 26 ) to give 



log(-/i 2 ) \og(K 2 K 2 /st) + dilog(l - g/jg) - dilog(l - Kl/t) 

st-K\K\ +UW - 



The expression (29) can be rewritten as 



r \ogs + \og(l + K 2 x/s) r°° \og(s + K 2 x)-\og(K 2 + tx) 

J (s + K 2 x)(K 2 + tx) X J (s + K 2 x)((s + K 2 x)-(K 2 + tx)) 

Of these terms, the first immediately yields 

\ogs\og(K 2 K 2 /st) + dilog(l - K 2 K 2 /st) 
st - K 2 K\ 
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and for the second term, perform the Mobius transformation 



w 



s + K 2 x 
Kl + tx 



(This choice turns on the fact that ((s + K 2 x) — [K\ + tx))' 1 is a removable 
singularity. It is tantamount to a proof of Abel's functional equation for the 
dilogarithm.) It then becomes 



st - K\K\ 



l-KSl* log(l - W) 



dw 



l-s/K* M 1 - W ) 

\ \og 2 {Kl/t) - \ \og 2 {s/Kl) + dilog(l - KUt) 



dilog(l - s/Kl) 



st - K\K\ 



(30) 



Combining all the terms, we recover the symmetry in K 2 and K\ in the for- 
mula: 



A 2me (s,t,KlKl^) = 

2 log(-// 2 /V^) \og{K 2 K 2 /st) + 2 dilog(l - K 2 2 Kl/st) \ 
-(dilog(l - K 2 /s) - dilog(l - 
-(dilog(l - tf 4 2 /t) - dilog(l - s/K 2 )) 
-(dilog(l - tfj/s) - dilog(l - t/i^|)) 
-(dilog(l - K 2 /t) - dilog(l - s/K 2 )) 



st 



KlK\ 



>(31) 



The complete symmetry is shown by the longer expression: 

\og^\stK 2 K 2 )-h)\og{K 2 K 2 /st) 
+ dilog(l - K 2 K 2 /st) - dilog(l - st/K 2 K 2 ) 
-(dilog(l - K 2 2 /s) - dilog(l - t/K 4 2 )) 
-(dilog(l - K 2 /t) - dilog(l - s/K 2 )) 
-(dilog(l - Kf/s) - dilog(l - t/K 2 )) 
-(dilog(l - K 2 /t) - dilog(l - s/X 2 )) 



st - K\K\ 



(32) 



As in the 3-mass integral, this formula shows how the case K\K^ = st is finite. 
Indeed, it may be seen that the integral allows both k\ and k 3 to be zero. In this 
extreme case there are only two distinct and only one kinematical invariant, 
s = t = K\ = K\. The integral is then log(/i 4 /s 2 )/s 2 . The method given above 
may be considered as using this case as baseline and seeing the effect of varying 
ki and k 3 away from zero. 
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7.3 O-mass case 



In the more degenerate cases, where adjacent edges of the Feynman tetrahedron 
lie on the singular quadric, the integral diverges as log 2 \x rather than as log /x, and 
more intricate work is necessary to separate the singular part. 

The most degenerate case is that of the O-mass integral. Then there are just four 
momentum twistors A, B,C, D. The parameters xi,X2,Xs,X4 correspond to lines 
AB, BC, CD, DA respectively. 



We may take from the preceding analysis the expression (25), and set K% = K\ = 
0. In this case there is no point in deferring the ^-integration, after which it 
becomes 

- 1 r y{2 - y) du 

o V Jo (us - 4ii?)((y(2 - y) - u)t - 4(1?) ' 
which again we split into two pieces again by writing 2 = 2y + (2 — 2y). 

The first piece is not now infra-red finite, but can still be easily evaluated as: 
2 f 1 2 log(2/(2 — y)) + log(- V/s ) + log(- V/t) dy + 



st Jo 2-y 

= 4«log2) 2 - tt 2 /12 - (log(-V/ S ) + log(-4/i 2 /t)) log 2 + 0(ji) . (33) 
st 

The second piece is: 

du dw 



st J Ja(u- 4h 2 /s)(w - 4/i 2 /f) 



4{log(-V/ S ) log(- V/t) - vr 2 /6} + GO*) . (34) 
st 



Combining the terms, we have 



A 0m (s, t, /i 2 ) = ^{2 log(-/i 2 / S ) log(-// 2 /t) - tt 2 } + 0(ai) . (35) 



7.4 1-mass case 

When K\ = 0, but K\ ^ 0, the first piece of the integral becomes 



rl r v(2-y) du 

2 / ydy 



'o 



o (y(2 - y)Kl + (t - K 2 2 )u - 4/i 2 )((y(2 - y) - u)s - 4/x 2 ) 
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Integrate out u, and obtain 

2_ f 1 dy 2\og(y(2-y))-\og(^Kl/st) | Q ^ 
s ~t Jo 2 — y 

= ^(log41og(-VK 2 2 / S t) + (log2) 2 -7r 2 /12) + 0(/i). (36) 
The second piece is 

, . du dw 

2 ' 



! A {us + wK\ - 4fi 2 ){wt - Aft 2 ) ' 
where A is the triangle < w,u,w + u < 1. Integrating out u, this is 
2 f 1 dm 

- J o ~ t _ 4 2 log(l - w(l - X 2 2 /s) - 4/i 2 / S ) - log(«,fff /a - 4/i 2 / S ) . 

The first logarithm is infra-red finite and integrates immediately to 

^dilog(l-A- 2 /s) + 0(/i). 

In the second logarithm substitute x = —Afx 2 /(wt — Afi 2 ) and hence obtain 
2 f 1 dr 

- / - {\og{- A^Kl/stx) + log(l - x(l - t/Kl))} + 0(jm) 

st J-AiP/t x 

2_ / (log(-4/i 2 /t) log(-4/x 2 ff*/st) - \\og\-^ 2 /t) 
st 1 +dilog(l-t/ J K- 2 )) + 0(/i) 



(37) 



These terms combine to give: 



21og(- / i 2 / S )log(- /U 2 A)-log^(-/i 2 /K 2 ) 
A lm (s,t,K 2 2 ^ 2 ) = -{ -2dilog(l-K 2 / S )-2dilog(l-K 2 /t)-7r 2 /3 >.(38) 

+0(/i) 



7.5 2-mass-hard case 



For the '2-mass-hard' case, with K 4 and K 3 non-null, we have six momentum 
twistors A, B, C±, C 2 , D 1 , D 2 . A little more ingenuity is required for this integral. 
The method given here is probably not optimal, but gives an indication of the 
geometrical relationships underlying the structure of the result. We start with 
the fully symmetrical Feynman tetrahedron as the parameter space, so that the 
integral is: 

f 1 f 1 f 1 f 1 + ^2 + «3 + «4 — 1) daida^dasda^ 

Jo Jo Jo Jo ("lM + a 2 a 4 s + a 2 a 3 K% + a 3 a 4 Kf - /i 2 ) 2 ' 
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Note that symmetry between x 2 and x 3 implies a symmetry between K 2 and K 2 . 
Eliminating and a±, this is 



'1 - a 2 - a 3 ) da 2 da 3 



(39) 



a ( (« 2 (1 — «2 — "3)5 + a 2 a 3 K 2 + a 3 (l — a 2 — "3)^4 - /-t 2 ) 
y (a 3 (l — a 2 — a 3 )t + a 2 a 3 K 2 - /^ 2 ) 

where the region of integration A is the triangle < a 2 , ot 3 < 1, a 2 + oc 3 < 1. 

Now this integral, but with iff = 0, has already been evaluated, as it is a 1-mass 
integral. We therefore study the difference between the 2-mass-hard integral and 
this 1-mass integral. This difference is 



a 3 K\ (1 — Qf 2 — « 3 ) 2 da 2 da 3 



a / (a 2 (l - a 2 — a 3 )s + a 2 a 3 K 2 - fx 2 ) 

(a 2 (l — oi 2 — a 3 )s + a 3 (l — a 2 — a 3 )if| + a 2 a 3 K^ - /i 2 ) 
(a 3 (l - a 2 — a 3 )i + c^a^Kf - /r 2 ) 



(40) 



Next, we use a splitting of the numerator factor 

cn 3 (l — a 2 — a 3 ) = [— a 2 a 3 K^/t] + [/i 2 /t] + [(a 3 (l — a 2 — a 3 )t + +a 2 a 3 K^ - ii 2 )/t] 

to divide this integral into three pieces. The first piece is infra-red finite, for up to 
0(/i) terms it is: 



-^3-^4 f f (1 — CH2 — a 3 ) da 2 da; 3 



A / ((1 - a 2 - a 3 )s + a 3 K 2 ) 

(a 2 (l - ol 2 - a 3 )s + a 3 (l - a 2 - a 3 )Kj + a 2 a 3 K. 
\ ((1 -a 2 - a 3 )t + a 2 Kl) 

To evaluate this, it is convenient to change co-ordinates to: 

y = (1 - a 2 - a 3 )/(a 2 K 2 ), v = (l-a 2 - a 3 )/(a 3 K 2 ), 

yielding 



(41) 



K 2 
ill 

t 



POO POO 

/ / (l + zt)- 1 (ys + K%z+l)- 1 (l + ys)- 1 dydz (< 
J Jo 

- / z-\l + zt)- 1 log(l + K\z) dz = --( dilog(l - K 2 Jt) - tt 2 /6) 
st J st 



by using equation ( 58 ) in the Appendix. 
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The second piece is 

fi 2 Kj (1 — a 2 — a 3 ) da 2 da 3 




(a 2 (l 


— «2 


- a 3 )s + a 2 a 3 K 2 - 


(a 2 (l 


— «2 


— a 3 )s + a 3 (l — a 2 


(a 3 (l 


- «2 


- a 3 )t + a 2 a 3 iff - 



(43) 



+0(Ai), (44) 



a {oi 2 s - /i 2 )(a 2 s + a 3 iff - fi 2 ){a 3 t - [i 2 ) 



which is essentially the same integral as (42), with /i 2 playing the role of iff, and 
again gives 



-i(dil0g(l-if 4 2 /t)-7T 2 /6). 



The third piece is 



if J f f (1 - a 2 — a 3 ) da 2 da 2 



t JJa( (a 2 (l - a 2 - a 3 )s + a 3 (l - a 2 - a 3 )K% + a 2 a 3 K% - {i 2 
\ (a 2 (l - a 2 - « 3 )s + ol 2 ol 3 K\ - /i 2 ) 

which we may write as 

iff ff , , /-i-«a-«3 dai 



(45) 



da2do; 3 



* 77a 7o (a 2 (l - a 2 - a 3 )s + aia 3 if| + 0205-^1 - /i 



2^2 



iff Z" 1 f 1 f 1 8(ai + a 2 + a 3 + a 4 — 1) d«i da 2 da 3 da 4 



t Jo Jo Jo Jo («2(«i + «4)s + a x a 3 K\ + a 2 a 3 iff - /i 2 ) 2 



= ^^(iff^iff,,. 2 ) = ^-A 2mh (Kl, s, iff, s, /i 2 ) (46) 

This object can be considered as the 2-mass-hard amplitude obtained from the 
one we are considering by merging ki and k 2 into a single non-null momentum 
ki + k 2 , but splitting if 4 into a sum of two null momenta, one of them orthogonal 
to ki + k 2 . 

In what follows it is convenient to use the dimensionless functions 

f 2 (s,t,KlKj) = stA 2mh (s,t,K 2 2 ,K 2 4 ), AM, iff) = stA lm (s, t, iff) , (47) 

which are just twice the conventional 'box functions'. In these terms (suppressing 
the dependence on \x and the 0{p) terms temporarily) we have the identity: 

f 2 (s, t, iff, if 2 ) = fi(s, t, iff) - f 2 (K I s, s, iff) - 2( dilog(l - if 2 /*) - tt 2 /6) (48) 
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By applying this identity three times we can find f 2 without any more integration. 
We have 



f 2 (Kis,s,K 2 3 ) 
UKlKlKls) 
f 2 (s,KlKlKl) 



h{Kl s, s) - f 2 (K 2 3 , Kl Kl s) - 2( dilog(l - K 2 Js) -tt 2 /6) , 
f x {Kl Kl K 2 4 )-f 2 (s, Kl Kl K 2 )-2{ dilog(l-s/^ 2 ) -tt 2 /6) 
his, Kl Kl) - f 2 (Kl s, s, Kl) - 2( dilog(l-K 4 2 /K 3 2 ) -tt 2 /6) . 



Hence 



2f 2 (Kls,s,Kl 



= h(Kl s, s) - 2( dilog(l - Kj/s) - vr 2 /6) . 

- h(Kl Kl Kl) + 2( dilog(l - s/Kl) - tt 2 /6) 

+ / 1 ( S ,^3 2 ,ir3 2 )-2(dilog(l-ir 2 /^3 2 )-7r 2 /6). (49) 



Using the known form of the f\ function, from (|38|), this yields 

2f 2 (Kls,s,Kl) = 



-4dilog(l - Kj/s) + log 2 (V/ S ) + log 2 (-/i 2 /K 2 
- log 2 (-/, 2 /Kl) - 2 log(K 2 3 /s) log(Kj/s) . 



(50) 



Hence 



A^s^KlKlv 1 ) 



1 

si 



{ 21og(- /U 2 / S )log(-/i 2 /t) -\\ g 2 (-^/s) 
-\\og 2 {-^/Ki)-\\og\-^/Ki) 
+ \og(Kj/s)\og{Kl/s) 
-2 dilog(l - K\lt) - 2 dilog(l - K 2 Jt) 

{ +0(/i) 



(51) 



Note that the symmetry in Kl and K\ is restored. It would be advantageous to see 
how these transformations and identities relate to geometric regions which break 
up the tetrahedron, and deal with the singular edge and vertex effects. 



8 Generalization to independent fii 

The 3-mass integral is the one degenerate case where it is simple to generalize 
from a single /i 2 to independent /i 2 , /u|, /X3, /i 2 in the obvious sense. We follow the 



same analysis as in section 7.1, and the only difference from (19) is that the term 



/i 2 (l + a + (3 + 7) 2 is replaced by (/if + a/1 2 + /3/i| + 7/14X1 + a + f3 + 7). 

By the same arguments the resulting integral may be simplified, up to terms which 
vanish as //j — > 0, to 

dad/3d7 

(a(K 2 + 7 s) + f3(t + 7 ^ 2 ) + a/3K 2 3 - (1 + 7 )(/i 2 + 7 /i 2 ))2 ' 
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Notably, /x 2 and pi 3 can be taken to zero; finiteness only requires regulators on 
those propagators which meet the null momentum. 



The integral will thus differ from (20) only in that log(/i 2 (l + r y) 2 ) is replaced by 



log((//i+7/i|)(l+7)), and so the result of the integral differs from the stated value 

d7 . 



of the 3-mass integral ( 22 ) by 

log((/i? + 7/il)/V(l+7)) 



(t + 1 Kl)(Kt + 1 s) 

This is a perfectly well-defined function of the kinematic scalars and /ii, /X2, but the 
generalization of our asymptotic formulas is more problematic, because statements 
about limiting functions depend strongly on exactly how the limits /ii, /14 — > are 
taken. If, for instance, the ratio /X4//X1 tends to some number A then this difference 
effect is 



1 



{log^W/x 2 ) \og(KlK 2 Jst) 

. - dilog(l - XKl/s) + dilog(l - Kt/s) 
st - A 4 K 2 ^ + dUog(1 _ Xt/K ^ _ dUog(1 _ t/K ^ 

We might like to make a comparison with the formalism of Alday, Henn, Plefka 
and Schuster (2009), which treats the effect of different mass parameters. But it 
is important to note that for different mass parameters their limit takes a quite 
different form from that discussed above. In their formalism, the external null 
momentum k\ is approached as a limit in which K\ = (//i — /U4) 2 . Such a limit is 
expressible in the formalism described in this note, but it is not very natural in a 
twistor-geometric setting, as there is no continuous transition from a small mass 
to zero mass. Either one momentum twistor is needed, or more than one. 



9 Summary 

The results given above are in agreement with standard results obtained by di- 
mensional regularization methods, as tabulated by Bern et al. (1993), Bern et al. 
(2004), in the following sense. First, allowance must be made for the fact that we 
have given the box integrals rather than the box functions, as described in section 
2 above. The correspondence may then be obtained by multiplying the dimension- 
ally regularized expression by (/i 2 ) e , subtracting simple and double poles in e, and 
then taking e = 0. 

We have neglected terms which vanish with /i 2 , but these have an important sin- 
gularity structure; they cannot be thought of as referring to a Taylor series in 
/1 2 . Their finiteness has the character of a finite value at a branch point, like 
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zlogz, log(l — z)\ogz or dilog(l — z) at z = 0. It would also be more correct 
to write them as 0(fi 2 /M), where M is the minimum of the kinematical vari- 
ables s,t,Kf, to indicate that the approximation of these divergent functions by 
functions involving logyU 2 is no approximation at all unless fi 2 is small compared 
with all the other parameters. This gives another way of seeing why the various 
limits /i 2 — >■ 0, K 2 — > do not commute, so that each degenerate case has to be 
considered separately. 

The various expressions for the box integrals show three levels of invariant struc- 
ture: there are divergent terms which are not even scale- invariant, there are terms 
like dilog(l — K 2 /s) which are scale-invariant but break conformal invariance, and 
there are conformally invariant terms like dilog(l — K\K\j st\ When the expres- 
sions are written in the natural momentum-twistor parameters, these features are 
immediately apparent. In the case of a six-field amplitude, for instance, when 
the expressions above are written in terms of six momentum- twistors Z\ . . . Z§, 
the term dilog(l - K 2 K 2 /st) becomes dilog(l - (1234) (4561)/(6134) (4512)) = 
dilog((1345)(1246)/(6134)(4512)), where (ijkl) = e^ZfZ^ Z%Zf , thus showing 
the manifest conformal symmetry. In contrast, the breaking of conformal invari- 
ance in the other terms is shown explicitly by their dependence on the infinity 
bi-twistor I. At this point we refer to the complementary work of Mason and 
Skinner (2010) which considerably develops this geometrical structure. Generally, 
it would appear that this exploration of the twistor geometry of the box functions 
gives only a first hint of new structures which are emerging from the analysis. 



10 Historical note 

Much interest in early twistor theory lay in the study of the integral 

(x - pi) 2 {x - p 2 ) 2 {x - p 3 ) 2 {x - p 4 ) 2 d X ^ 

arising in first-order 4 theory, i.e. in x-space, not in momentum space. Here the 
x-integral is over real x, and the Pi are four points in CM, of which two must 
be in the future tube and two in the past tube for a non-singular, non-vanishing 
result. Thus, l/(x — pi) 2 is a classical positive or negative frequency solution of 
the massless scalar field, called an 'elementary state' by twistor theorists. The 
choices between future and past tube gives rise to three 'channels', which cor- 



respond to the three logarithmic period functions noted in (17). (For instance, 
the first of those logarithmic functions allows a = c, b = 0, which corresponds to 
allowing pi = p 2) pz = Pa-) The first twistor-integration results were stated by 
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Penrose and McCallum (1972); curiously, these mis-identified the analogue of the 
'leading singularity' contour as the correct ('logarithmic') contour for an ampli- 
tude. Thus even in the earliest days, these questions of connecting contour with 
physical amplitude played a dominating role in twistor integration. This was soon 
corrected, as later surveys, such as by Hodges and Huggett (1980), Hodges (1983), 
indicate. 

The construction of the dilogarithmic result for the 4-mass integral by an S 7 x S 1 
contour appeared in an informally published note (Hodges 1977). Its motivation 
was different from that presented in this note; its title was 'Crossing and twistor 



diagrams' and its hope was to define a 'super-amplitude' for the 4 integral (52), 
in such a way that each channel would arise as a period. The construction used 
the linearity of CP 5 , leading to the basic dilogarithmic answer, and its relationship 
to the logarithmic periods. The statement of the result did not fully address the 
problem that the dilogarithmic result can only be regarded as a limiting value as 
elements of CP 5 approach the singular quadric. The original 'crossing' idea could 
not be generalized from this application to elementary states, and was never taken 
any further. This present note arose in the context of the renaissance of twistor 
field theory in the wake of Witten's twistor string model in 2003. It then became 
notable that (a) the dilogarithmic form of the 4-mass box integral is just that of 
the structure found in 1977 and (b) that the most obvious way to make the 1977 
construction rigorous is to add the fi 2 I al3 terms, and regard the finite result as a 
limit as fi 2 — > 0. 

Use of a generalized 'Feynman trick' was widespread in early twistor theory work, 
following Penrose's use of it in many examples. As an example, the study of the 
4 scalar integral (52) in (Hodges 1983) used the fact that the product of fields 
l/(x — pi) 2 (x — P2) 2 can be written as a Dirichlet average of fields of the simpler 
form l/((x — z) 2 ) 2 , where z varies over a sphere in CM. 

The 1977 note also used the fact that the contour for the fundamental twistor 
integral ^ need not be realised as an S 7 x S 1 . Instead, it integrated W a as a 
S 1 overall phase, times a tetrahedron in projective W a space, with vertices at 
W a = A a , B a ,C a , D a . The result is 

{27ii)~ i je a ^ 5 A a BpC 1 D & {{Z a A a ){Z a Bp){Z^C^{Z 5 D 5 ))- 1 d A Z. 

A direct product of four S 1 integrals in Z completes the integration. This contour 
is actually homologous to the S 7 x S 1 when four extra regions are added, each 
of which is inside a space of form W a = constant, and so contributes nothing to 
the integral. The use of this contour makes the basic Feynman loop integral ([7]) 
trivial, if the tetrahedron is chosen with vertices on O and /. 
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The motivation for introducing this complication was that whilst the S 7 x S* 1 con- 
tour clearly did not exist in the singular limit, some adaptation of this boundary- 
defined contour might in fact survive in the limiting case. This idea may not 
be completely wrong, in view of the central role of boundary-defined integrals in 
momentum- twistor space (Hodges 2009). In any case, it is notable that contour 
integration problems of the 1970s are very relevant today. 
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12 Appendix 



The dilogarithm function dilogz, also written Li 2 (z), is defined by 

ddo^-fi^d^f^. ,53) 



By elementary integration, 



dilog(l - = -dilog(f-^)-^log 2 z, (54) 



2 



7T 2 



dilog(l — z) = — dilogz + — — log(f — z)\ogz (55) 

6 

7T 2 f 

di\og(z l ) = -dilog^; - — - -\og 2 (-z) , (56) 

o 2 

showing how dilogz behaves under cross-ratio transformations. At z — 1, dilogz 
is not analytic, and has a branch point. But it also has the finite value 7r 2 /6 there. 
It is useful to consider dilog z as defined on the complex plane with a cut from 1 
to oo along the real axis; it then has a discontinuity of 2rri logz across this cut. It 
follows that the definite integral 



Jo 



log(l + Az) ^ 



! (z + k)(z + k) 

is equivalent to the contour integral 

1 £ dilog(l + \z) dilog(l — \k) — dilog(l — \R) 

~dz = z 7 i J 



2%i J (z + k){z + k) k — k 

and we use this repeatedly. Here k and k are assumed to be away from the cut. 
But by taking the limit k — > 0, we deduce 



Jo 



lQg(l + \Z) d ^ _ dilQg(l - \K) - 7T 2 /6 

(z + k)z k 
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